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Abstract 

The aim of this work is to formulate for dissipative system a least 
action principle that can keep all the main features of the conven- 
tional Hamiltonian/Lagrangian mechanics such as the Hamiltonian, 
Lagrangian and Hamilton- Jacobi equations, three formulations of the 
classical mechanics, with a unique single Lagrangian function which 
has the usual close energy connection and Legendre transformation. 
We consider a whole isolated conservative system containing a damped 
body and its environment, coupled to each other by friction. The La- 
grangian is L = K — V — Ed with an effective conservative Hamiltonian 
H = K + V + Ed where K is kinetic energy of the damped body, V 
its potential energy and Ed the negative work of the friction force. 
This Least action principle can also be derived from the virtual work 
principle. It is shown that, within this formulation, the least action 
principle can be equivalent to a least dissipation principle for the case 
of Stokes damping or, more generally, for overdamped motion with 
more kinds of damping. This extension opens a possible way to study 
the general relation between two independent classes of variational 
principles, i.e., those of Hamiltonian/Lagrangian mechanics and those 
associated with energy dissipation. 

Keywords: Classical mechanics, Variational principle, Lagrangian, Ac- 
tion, dissipative systems 
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1 Introduction 



The Least Action Principle (LAP) is one of the most valuable heritages 
from the classical mechanics [1, 2, 3]. The fact that the formulation of the 
whole classical physics as well as of the quantum theory in its path integral 
formalism [ 1] could be based on or related to this single mathematical rule 
gives to LAP a fundamental priority to all other visibly different principles, 
empirical laws and differential equations in different branches of physics. This 
priority of LAP has nourished two major hopes or ambitions of physicists. 
The first one is the (rather controversial) effort to deepen the understanding 
of nature through this principle and to search for the fundamental meaning 
of its exceptional universality in physics[3, 5, 6, 7]. The second one is to 
extend it to more domains such as thermodynamics and statistical mechanics 
(with the pioneer effort of Boltzmann, Helmholtz and Hertz [8]), stochastic 
dynamics (e.g., large deviation theory[9] and stochastic mechanics[10]), and 
dissipative mechanical systems[ll, 12, 13]. The present work is following this 
last effort to formulate LAP for dissipative or nonconservative mechanical 
systems. 

LAP was originally formulated only for Hamiltonian system[2], i.e., the 
sum H = K + V of kinetic energy K and potential energy V of the consid- 
ered system satisfies the Hamiltonian equations. For Hamiltonian systems, 
any real trajectory between two given configuration points must satisfy the 
LAP given by a vanishing first variation 5 A due to tiny deformation of the 
trajectory [2, 6]: 

SA = S I '" Ldt = f" SLdt = (1) 
Jo Jo 

where the action A = J b Ldt is a time integral of the Lagrangian L = K — V 
on the trajectory over a fixed time period tb- One of the important results of 
this variational calculus is the Euler-Lagrange equation given by[6] (for one 
freedom x) 

d (dL\ dL 
dt \dx J dx 

where x is the velocity. In many cases when H and L do not depend on time 
explicitly, a Hamiltonian system is energy conservative. The problem comes 
from nonconservative systems having, say, a dissipative force fa which must 
be introduced into the equation of motion in this way ^ f ||J — = fd- This 
is equivalent to J tb (SL + f d 5x)dt = 0[1 1]. Despite this vanishing equality, it 
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is impossible to calculate and optimize an action integral like A above with 
a single (Lagrangian) function satisfying Eq.(2). This difficulty leads to the 
disappearance of LAP in dissipative systems. 

There has been a longstanding effort to formulate LAP for nonconserva- 
tive or dissipative system[12]. As far as we know, the first proposition was 
made by Rayleigh[16] who introduced a 'dissipative function' D = \Cx 2 to 

write ^ {jji) + §f — = 0' where C is the drag constant in the Stokes law 

fd = —m^x and m the mass of the damped body. Although the equation of 
motion is kept in a similar form as Lagrangian equation, LAP is not recov- 
ered since there is no single Lagrangian for defining an action which satisfies 
Eq.(f). Other major propositions include the Bateman approach [17] to in- 
troduce complementary variables and equations, the definition of dissipative 
Lagrangian by multiplying the non dissipative one with an exponential factor 
exp((t)[lS] where t is the time, the fractional derivative formulation[19], and 
the pseudo-Hamiltonian mechanics [20] where a parameter was introduced to 
characterize the degree of dissipation. The reader is referred to the reviews 
in [11, 12, 15, 19, 20] about the details of these propositions. In general, the 
Lagrangian in these solutions is not unique and has no energy connection 
like L = K — V (see for instance the quasi-Lagrangian L = e&(K — V) and 
the corresponding quasi-Hamiltonian H = e~&K + e&V for damped har- 
monic oscillator [18]). Hence no variational or optimal calculus was possible 
in general formfll, 12, 15]. 

A common character of these works is that the damped body is the only 
object taken into account in the calculations as if it was isolated. However, a 
dissipative system is always coupled to an environment and loses energy into 
the latter, an integral part of the motion. As far as this lost energy is not 
considered, the quasi-Lagrangian function of the damped body inevitably 
loses energy connection and generic optimal charactersfll, 12] as mentioned 
above. 

The aim of this work is to establish a LAP of dissipative system that 
recovers the energy connection and the uniqueness of a single Lagrangian 
function, its relation with a conservative Hamiltonian (Legendre transform), 
as well as the three formulations of analytical mechanics, i.e., the Hamilto- 
nian, Lagrangian and the Hamilton- Jacobi equations. 
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2 The conservative Hamiltonian 



The idea is to consider the damped moving body and its environment, cou- 
pled to each other by dissipative force, as a whole conservative system to 
which LAP can be a priori applied. The total Hamiltonian includes the 
instantaneous kinetic and the potential energy of the body, as well as the 
mechanical energy that is lost from the beginning of the motion and trans- 
formed into heat or other forms of kinetic and potential energy (noises, vi- 
bration etc) in the environment. Concretely, we construct a total system 
composed of a one dimensional large moving body (system 1) along the axis 
x and its environment (system 2) which includes all the parts coupled to 
system 1 by friction and receiving the dissipated mechanical energy. The en- 
ergy transfer from system 1 to system 2 occurs only through a friction force 
given, for example, by the Stokes law. The total Hamiltonian can be given 
by H = K\ + Vi + H 2 where K\ = ^mx 2 is the kinetic energy, V\ the potential 
of a conservative force acting only on system 1, and H 2 the energy of system 
2. The interaction energy between 1 and 2 is neglected or supposed to be 
constant and dropped from the variational calculus. For simplicity, suppose 
that system 2 does not move as a whole, hence the macroscopic moving paths 
of the whole system are just the paths of system 1. This allows to calculate 
the action of the whole system along the paths s = x(t) of system 1 moving 
between two configuration points x a and Xb during the time period % — t a 
(t a = from now on). The amount of energy Ed dissipated from system 1 to 
2 from t a = to a time moment t (0 < t < t^) is given by the negative work 
of the friction force: 

rx(t) _ fx(t) ft 

Ed[x(0,t)} = — / fd ■ ds(r) = m( / x{r)dx{r) = m( / x 2 {r)dT (3) 

JXa JO JO 

where r is any time moment between t a = and t, ds = dx = xdr a small 
displacement along s. E d [x(0,t)} means that the dissipated energy depends 
on both the past trajectory from x a = x(0) to the instantaneous position 
x(t). 

On the other hand, during the motion, the energy of system 2 at time 
t can be written as H 2 = H 2 + Ed where H 2 is its energy (a constant) at 
t a and can be dropped from the variational calculus of the action. Finally 
we will consider only the following effective Hamiltonian for the motion : 
H — K i + Vi + Ed- This effective Hamiltonian is formally nonlocal due to 
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the space-time non locality of the integral of E d in Eq.(3). E d is dissipated 
part of H 2 , its non- locality comes from its expression on the coordinates x 
of the damped body. But H2 is actually a local function of the motion and 
can be expressed by the instantaneous energy of the N constituent particles 
of the environment at the moment t, i.e., if 2 = J2f=i(ki + Vi) where ki and Vi 
is respectively the instantaneous kinetic and potential energy of the particle 
i. 

Anyway, this non-locality of the Hamiltonian can have some influence 
on the variational calculus. It will be shown later that this influence can 
be avoided by the consideration of the principle of locality or of the energy 
conservation. The instantaneous increment of E d is compensated by the 
simultaneous equal decrement of K1 + V1, assuring a constant H in time and 
space for the isolated whole system. 

Before proceeding with Lagrangian function and LAP, we stress that the 
impact of the thermal fluctuation in system 2 on system 1 should be ne- 
glected in order to have a smooth and deterministic motion of the latter. 
This is not difficult for a body which is much larger than the constituents of 
system 2 and has much larger energy variation during the motion than the 
energy fluctuation of the thermal motion in system 2. The reason for this 
approximation is that LAP for stochastic motion is still an unsolved problem 
to date. LAP in its conventional form is not compatible with random dy- 
namics. It is also for this reason that, in this work aiming at extending LAP, 
expressing E d in terms of the coordinates and velocities of the microscopic 
constituents of system 2 will not help because of their random motion re- 
lated to heat, an inevitable effect of friction. Fortunately, no matter how the 
particles in system 2 move and whatever is their state of randomness, their 
unique relevant effect on the motion of system 1 is the friction. Thanks to 
the conservation law of energy, the dissipated energy E d , i.e., the increment 
of the sum of the kinetic and potential energies of all the particles in system 

2 expressed in terms of their coordinates and velocities, whatever is the form, 
can be mapped onto the coordinates of the system 1 through the work of the 
friction forces. 

3 A dissipative Lagrangian function 

The first difficulty for writing the Lagrangian is the impossibility to separate 
E d into kinetic and potential parts relative to the coordinate x of interest. 
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The second difficulty can be explained as follows. Suppose the simple case 
where Ed is only kinetic energy (heat in an ideal gas for example), we are 
inclined to write L = Ki—Vi+Ed according to the convention. Unfortunately, 
it is straightforward to show that this will lead to an incorrect equation of 
motion when this L is introduced into Eq.(2). 

In what follows, we consider the fact that Ed, after the integration in 
Eq.(3) over a certain trajectory from x(0) to the instantaneous position x(t) 
which is changing in time, implying that the integral of Eq.(3) is an indefi- 
nite one. According to the second fundamental theorem of calculus[21], the 
friction force at time t can be calculated from Ed by fd = — ^fjr in a similar 
way as the conservative force is derived from a potential. Obviously Ed is not 
a potential since it depends not only on x(t), but also on the past trajectory 
along which the integral Eq.(3) has been carried out (in practice, a trajectory 
s — s(t) can be introduced in the calculation of Ed by writing ds = s(r)dr). 
Moreover, it is impossible, contrary to potential energy, to recover Ed (an 
increasing function of time) as mechanical energy of system 1 just by moving 
the latter backwards. However, Ed has an common character with potential 
energy : its instantaneous increase yields the resistance force through the 
instantaneous displacement of the body; in other words, the infinitesimal in- 
crease dEd at time t is equal to the negative work done by the friction force 
over an instantaneous displacement dx(t), i.e., dEd — —fd • dx = —fddx (for 
more than one dimension, this means dEd = —fd • dr or fd = — Vf^-E^)- We 
think that the above arguments are sufficient, from the energetic point of 
view, for considering Ed as a pseudo-potential and writing L = K\ — V\ — Ed 
as an effective Lagrangian. The effective action of the whole system on a 
given path between a and b is then given by 



A = / (K x -Vi- E d )dt = / (-mi 2 -Vi-mC x{r)dx(T))dt. (4) 



4 Variational formulation 
4.1 A global variational calculus 

Notice that this Lagrangian L is nonlocal since it depends on the state at 
time t as well as the state at time r, making A a non local action. This may, 
at first glance, complicates the variational calculus since, mathematically, 
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Figure 1: Illustration of an variational path (dotted line) over the entire 
optimal path (thick line) from the point a to the end point b. 



the variation of the action of Eq.(4) due to the variation of the whole path, 
illustrated in Fig.l, depends on both 8x(t) and Sx(r) at two different time 
moments. However, from physical point view, the consideration of the contri- 
bution of 5x(t) to 5 A is equivalent to considering the effect of the dissipative 
energy of the anterior motion on the present motion. This would violates 
the principle of locality of classical mechanics. The fact is that this energy 
has been already absorbed irreversibly into the surroundings and should not 
affect the mechanical motion any more if the thermal fluctuation of the en- 
vironment is neglected as mentioned above. Hence we propose to consider 
only the contribution of 5x(t). The deformation of the whole path is taken 
into account through the time integral over t from t = to t = T. The usual 
variational calculus [2, 6] is 



5A= [ SLdt 



o 



dL . , dL . . . , 
-5x{t) + —-^Sx(t) 



dx(t) dx(t) 



dt, (5) 



which becomes 8A = Joi'^m ~ iKdJ^))^ 2 ^)^ after the time integral by 
parts of 5x with the condition Sx(t a ) = 5x(tb) = 0. Since 5x(t) is arbitrary, 
the vanishing first variation 5 A = necessarily yields the Euler-Lagrangian 
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equation Eq.(2) and the Newtonian equation of damped motion: 



mx 



<W + E d ) 



dVi 



(6) 



dx dx 

where 'x is the acceleration. This is the Newtonian equation of system 1 
with Stokes friction. But here it was derived by vanishing variation of the 
action of the whole system, meaning that it should be considered as an 
equation of motion of the system 1+system 2 mapped on the coordinates of 
the system 1. Notice that, although the Stokes law was used in the above 
calculations, the result should be valid for any friction force f d = f(x) = 
-£lo {t) mr))dx(r). 



4.2 A local variational calculus 

In the above calculus, we considered the principle of locality of classical 
mechanics in order to avoid the influence of earlier states on the present 
motion. This help from another principle to the variational calculus is in fact 
not necessary if we consider the differential version of LAP. The argument is 
the following. If A is a minimum over the entire optimal trajectory between 
a and b, the same must be true over any segment of the trajectory, i.e., the 
time integral of L over a segment Ax must be a minimum whatever its length 
is. If not, we can always play with this segment to make A smaller than its 
minimal value along the optimal path. Now if AA is the action over this 
small segment around the time moment t, from Eq.(4), we have 



AA 



t+At I 



x(t) 



[-mx 2 — V\ — m( / x(r)dx(r))dt. (7) 
n 2 Jo 

The variation of the first two terms in the integrand is the same as in Eq.(5). 

The variation of the third term, i.e., of E d , is 5E d = ( j®w +Sx w x(r)dx(r) — 

m( Jq x(r)dx(r) = -^f^5x(t). Put this back into the variation of AA, we 

get 



5AA 



r t+At d 

t dt 



dL 



■8x(t) 



dt + 



r t+At gi 



d , dL 



dx{t) 
dL(t + At)^ . 

rt+^K dL d , dL 



8L(t) 



[ dx{t) dVdx{t) 



dx(t) 
)]8x(t)dt 



[ dx{t) dVdx{t) 



5x(t) 



)}5x(t)dt(8) 



(9) 
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Since the variational path must join the optimal path before and after the 
small Ax(t), hence 5x(t + At) = 5x(t) = 0. The Euler-Lagrangian equa- 
tion Eq.(2) will be a necessary consequence of the vanishing first variation 
5AA = 0. This confirms the LAP by the variation of whole path aided by 
the principle of locality. 



4.3 Derivation from virtual work principle 

The above LAP can also be derived from other fundamental principles as 
has been done in analytical mechanics by using virtual work principle of 
d'Alembert[2, 22]. This latter principle is valid in the presence of friction 
force. For 1-dimensional moving body, it reads: 

5W = (/i +f d - mx)5x = (10) 

where f\ = — ^ is the conservative force. This expression implies that 

system 2 has been involved in the motion since E d = Jq® f d dx(r) is a part of 
its energy. Using xSx = j- t (x8x) — xSx = ^(x5x) — ^±5x, then integrating 
Eq.(10) over time from t a = to t b , we get 

/ (-^5x -±5x)dt= (—Sx+—8x)dt= 5Ldt = 5A = 

Jo ox ox Jo ox Ox Jo 

(11) 

where we used 5x(t a ) = 5x(tb) — 0, L — K\ — V\ — E d and A = J <b Ldt. 

Within this formalism, it is easy to verify that the Legendre transforma- 
tion H = px — L is still valid, where p is the momentum of system 1. With 
the usual method[6] using Euler-Lagrange equation Eq.(3), the Hamiltonian 
equations can be derived : 

p=- d —x= d — (12) 
dx ' dp 

The Hamilton- Jacobi equation also holds. To see this, we relax tf, in the 
integral Eq.(4) or consider the integral as indefinite, and compute L = ^ = 
c + 7j^. Thanks to Eq.(2) and the above Legendre transformation, we can 



dA 

dx 1 dt 



get P = and the Hamilton- Jacobi equation : 

dA , . 

or ^ = — ^(ff ) 2 — V\ — E d to show the dissipative character of this equation 
for the whole system (damped body + environment). 
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4.4 Application of Maupertuis' principle 



It is worth mentioning that, if the conservation of the total Hamiltonian H 
is considered as a constraint of variation, the mathematical trouble of whole 
variation with the non locality of the expression of Ed can be easily avoided. 
The constraint of constant H for the total isolated system is reasonable be- 
cause any conceivable motion (even virtual) should not violate this universal 
law. With this in mind, the non locality of Lagrangian disappears if we con- 
sider the Legendre transformation L = px — H since L varies in the same 
way as px if H is constant. Hence the optimization of L is equivalent to that 
of the function px. This remind us of the Maupertuis least action principle 
with the action defined by A M = f* b mxdx[l}. It is known[15] that the LAP 
5 A = is equivalent to the Maupertuis principle stipulating 5A M = with 
the conditions of fixed points a and b as well as fixed total energy H instead 
of tb. The calculus is 

rb r b 



no no 

5A m = 5 / Pdx = / (5pdx + p5dx). 

J a J a 



Considering dx = xdt in the two terms, the first term becomes 8{^—)dt. Then 
making integration by part of 5x in the second term, considering 5x(a) = 
5x{b) = 0, this term becomes —mxbxdt. The total energy conservation 
8H = means 5(£—) = — ^5x — ^M^Sx. Finally, we have 

2m I ox ox J 1 



5A 



m 



dV 1 dE d 

mx 



dx dx 



8 xdt, 



which implies that the Maupertuis principle 5A m = necessarily leads to the 
Euler-Lagrange equation Eq.(2) and to Eq.(6) as well. 



4.5 Optimization of dissipation? 

A notable aspect of this dissipative LAP is the inclusion of dissipative energy 
in the variational optimization through the relation between L (or action A) 
and Ed (or its time integral Ad = Jg 6 E d dt). This naturally makes us think 
about possible relationship between this LAP and the optimization of en- 
ergy dissipation. The conventional LAP and the variational principles of 
energy dissipation are two independent families of axioms (there are even 
assumptions combining LAP and a dissipation principle as two independent 
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principles which hold simultaneously for a dissipative fluid system [23]), each 
being valid for its own systems. Many variational principles of dissipation 
have been formulated in relation with thermodynamic phenomena and en- 
tropy production [11, 24]. Few attention has been paid to extremum prin- 
ciples of dissipation in mechanical motion. An example of this is the path 
of least dissipation which, rather by intuition, states that a mechanical sys- 
tem should follow the path of least energy dissipation in the case of damped 
motion. But in view of the present LAP, this principle is not that evident. 
From the vanishing variation S f b (Ki — V± — Ed)dt = 0, it is obvious that 
8 Ad 7^ 0. Hence in general there is no extremum or stationary of Ad on 
the path of least action if we do not change the conditions of the vanishing 
variation 5 A = 0, i.e., fixed time duration tf, and fixed initial and final points 
a and b. These conditions can be of course modified for different variational 
problems. It is instructive to see the case of Maupertuis action. It is easy to 
verify that, in the case of Stokes friction, the Maupertuis action can be writ- 
ten as Am = ^ m(xdx = where E d = f** (xdx is the negative work of 
the friction force over the entire trajectory from a to b. Hence Maupertuis 
principle implies SE d = 0. In other words, the path of least action is just the 
path of least resistance in the case of Stokes friction. 

Finally we would like to mention that, although the term "least action" 
is used here for historical reason, the stationarity 5 A = is not necessarily a 
minimum. The nature of the stationarity (minimum, maximum or inflection) 
of A has been addressed in our recent work [30] by numerical calculation of 
the action A op along the optimal path compared to the actions of numer- 
ous other paths created with tiny deformations of the optimal one. The 
damped body is a small ball subject to different conservative (constant and 
harmonic) forces combined with different (constant, Stokes and quadratic) 
drags, respectively. Fig. 2 gives an example of the result for a motion damped 
by Stokes friction in the gravitation field on the earth. The ball moves during 
1 s. The paths are discretized into n s = 1000 steps of equal time increment 
in the simulation. The reader is referred to [30] for further details. The 
conclusion is that the stationarity of A op undergoes a transition from mini- 
mum to maximum when ( and the dissipated energy increase from zero to 
very large value (overdamping) . This maximum dissipative action may have 
a hint about the optimization of energy dissipation. Let us make a varia- 
tion of the optimal path by a tiny but very intensive zigzag motion close 
to the optimal path. In the expression A = ^{Ki — V\ — EAdt, the fu- 
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tegral Jq 6 V\dt almost does not change because the variational path is very 
close to the optimal one, but the integral A k = Jq 6 K\dt and A d J^ b E^dt will 
increase enormously due to the larger variational velocity caused by the vio- 
lent zigzag motion. Hence the maximum action implies A over the variational 
path must be smaller than A op . In other words, Ad must increase, and more 
quickly than A^. This is a minimum energy dissipation over the optimal path 
of overdamped motion. Further careful work is necessary to elucidate other 
cases for different conditions of variation and different extremum principle of 
energy dissipation [25, 26]. 

5 Concluding remarks 

In summary, we formulated a possible answer to a longstanding question 
of classical mechanics about the least action principle for damped motion, 
in keeping all the four conventional formulations of mechanics, i.e., Newto- 
nian, Lagrangian, Hamiltonian and Hamilton- Jacobi equations. This work is 
based on the model of a conservative system composed of the moving body 
and its environment coupled by friction. It was shown that this system with 
"internal dissipation" satisfies both Lagrangian and Hamiltonian mechan- 
ics, leading to correct equation of damped motion in a general way. It was 
also shown that, within this formulation, the Maupertuis principle is equiv- 
alent to a least dissipation principle in the case of Stokes damping. A more 
general least dissipation principle is also discussed for the overdamped mo- 
tion. We hope that these results are helpful for further study of the relations 
between the variational principles of energy dissipation and the fundamen- 
tal principles of Lagrangian and Hamiltonian mechanics (see the efforts for 
stochastic dissipative systems in, for example [27, 28, 29]). It is also hoped 
that the present result is useful for the study of quantum dissipation in view 
of the role of action in the quantum wave propagator ip = e tA ^ h [l] and the 
close relationship between the Schrodinger equation and the Hamilton- Jacobi 
equation[16]. 
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Number of paths 
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Figure 2: Illustration of the transition of extrema by comparison of the action 
of the optimal path (dots) with the actions of other paths (circles) created 
by tiny deformations of the optimal one. (a) For £ = 0.1 s" 1 or smaller, A 
is in the bundle of smallest actions, (b) For £ = 1 s _1 , A, 



*-op 



1 or smaller, ^ op 
is in the middle 



rank, (c) For £ = 10 s 1 or larger, A op is in the bundle of largest actions. 



(d) ( dependence of the quantity AA 



A-A 



where A is the average 



\A \-\-\A p | 

action over all the paths. AA can be used to characterize the evolution of 
the extrema of A in three regimes: the minimum regime (AA > 0) for small 
C, the maximum regime (AA < 0) for large £ and the saddle point regime 
around A A = corresponding to a critical £ c = 1. The steep increases at the 
two extremities of the ranking are due to the insufficient number of paths 
around the smallest and largest actions. 
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